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Nonlinear evolution equation at small x with impact parameter dependence is analyzed numeri- 
cally. Saturation scales and the radius of expansion in impact parameter are extracted as functions 
of rapidity. Running coupling is included in this evolution, and it is found that the solution is sen- 
sitive to the infrared regularization. Kinematical effects beyond leading logarithmic approximation 
are taken partially into account by modifying the kernel which includes the rapidity dependent cuts. 
While the local nonlinear evolution is not very sensitive to these effects, the kinematical constraints 
cannot be neglected in the evolution with impact parameter. 



I. INTRODUCTION 



o 
O 

\l I High energy limit in Quantum Chromodynamics is one of the most intriguing problems in strong interaction physics. 

Since the Large Hadron Collider has recently opened a new kinematic regime, it is of vital importance to understand 
(~| ,; how one can calculate the cross section at high energies from basic principles in QCD. The BFKL equation obtained 
in the Regge limit of QCD [ll-13| predicted fast growth of the cross section with the energy, due to the exchange 
^ of the hard Pomeron. This behavior is however known to violate the unitarity bound of the scattering amplitude. 
^ Higher order corrections to the BFKL tame the growth but other terms due to the parton recombination 

, and multiple Pomeron exchanges have to be taken into account in order to guarantee the unitarity of the scattering 
amplitudes. 

^ ' The Balitsky-Kovchegov equation was derived in [l6j in the dipole approach \yf\ to high energy scattering in 
^ QCD and independently in the formalism for the operator product expansion at high energy as an evolution of the 
^-H ' Wilson line correlators with respect to the rapidity |18l - l22| . This equation also emerges in the Color Glass Condensate 
: formalism as a limit of the JIMWLK functional equation (23l - [27| , see also [28l. . The BK equation was shown to 
be equivalent to the BFKL evolution with the triple Pomeron vertex [sol . l3ll| when the latter one is restricted to the 
Mobius space of functions ^32]. 
• Numerous analyses of BK equation which were performed up to date, see for example [33l - l37j . focused on finding 
^-H I the solutions to this equation under the assumption that the amplitude does not depend on the impact parameter but 
only on a dipole size and rapidity. This assumption makes the equation relatively easy to solve (at least numerically). 
' ^ [ Usually one is justifying this approximation by considering the scattering off a very large target, therefore bringing a 
^ , large scale into the problem. This results in breaking the symmetry between the infrared and ultraviolet regions even in 
the case of the fixed coupling by generating a saturation scale which acts as a boundary and suppresses the diffusion 
into the infrared region ^35i . i38j . In the momentum space this approximation is equivalent to taking the forward 
' limit in the evolution with the simplifying approximation on the triple Pomeron vertex resulting in zero momentum 
transfer through the vertex in both reggeized gluon pairs [39|. On the other hand, the kernel of the equation has the 
property that it is invariant with respect to the translations, dilatations, rotations and inversions, which are Mobius 
transformations. Therefore the solutions to the nonlinear equation should also reflect this symmetry properties, at 
least in the leading-logarithmic, fixed coupling limit. On a deeper level it is related to the Mobius invariance of the 
2—^4 reggeized gluon transition vertex which was demonstrated explicitly [l^l. In the relation to the experiment, 
the detailed knowledge about the dynamics and expansion in impact parameter is of the utmost importance for the 
qualitative and quantitative description of the multiparticle production in hadronic collisions. 

In this paper we analyze this equation numerically including the full impact parameter dependence. This is an 
extension of the previous work [41] where this type of analysis was performed, see also [l^]. We significantly improve 
over [4l| the numerical accuracy and technique which enables us to evolve the equation much faster and to very 
high values of rapidity, of the order of ~ 50. In this way one can more accurately extract the asymptotic values 



•Electronic address; jxbl024@psu.edu] 

t Electronic address: ast asto® phys . psu. edu | 



2 



of different exponents which govern the growth of the saturation scales in this equation. We confirm the results of 
[4l| . on the dependence of the scattering amplitude as a function of dipole size and demonstrate that it vanishes 
for large dipole sizes. We also find the fast diffusion of the solution in impact parameter space and recover the 
power tails. The saturation scale both for small and large dipoles is extracted, and the dependencies on the impact 
parameter and rapidity are found. The results of the solutions to the equation in the leading logarithmic approximation 
(LL) are compared with the modified version of the equation proposed in 43]. The modified version contains the 
cutoffs in rapidity which originate from kinematical constraints. These cutoffs contain kinematical constraints in 
only approximate way but we know from the analysis of forward BFKL in momentum space that these constraints 
are known to reduce the speed of the evolution in a significant w av H . (for a related analysis on impact parameter 
dependence in nonlinear equation and the energy conservation see [45| . The BK without impact parameter dependence 
and with rapidity cutoffs was also analyzed in [ssl . |46| ) . We also include running coupling in our analysis and find 
that the effect of the running coupling is quite different than in the case without the impact parameter. In this paper 
we consider a prescription for the running coupling with the external dipole as the scale as well as the prescription 
derived in [l^]. The impact parameter dependent equation is extremely sensitive to the large dipole sizes and this is 
the region where the running coupling is very large and needs to be regularized by some other mechanism. 

In this analysis we did not attempt to regularize the large dipole size region in any way. It is at present totally unclear 
how confinement effects should be consistently included in the dipole formalism. Of course, for any phenomenological 
applications such cut should be included, perhaps similarly to what was done in [iSij . As we were interested in general 
properties of the evolution we did not attempt here to introduce additional cuts on large dipole sizes (via masses), 
which would interfere with the specific dynamics of the evolution. 

The paper is organized in the following way. In the next section, Sec. |TT] we briefly present the BK equation and 
discuss the modified version which includes the cutoffs in rapidity. In Sec. IIIII we describe the numerical methods of 
finding the solution. In Sec. II VI we first show the results of the solution without the impact parameter and extract the 
saturation scale for both the LL and the modified equation. In Sec.|V]we present the solutions with impact parameter. 
We discuss the form of the amplitude as a function of the dipole size, extract the saturation scales (both for small and 
large dipoles), and discuss the form of the impact parameter profile which emerges in the evolution. We present the 
solutions both in the case of the LL and for the modified kernel. Using the representation in terms of the conformal 
eigenfunctions we discuss the origins of different peaks in the amplitude as well as present estimates for the rapidity 
dependence of the small and large dipole saturation scales and the expansion radius in impact parameter. We also 
present the estimate of the cross section of the black disc radius and its dependence on the rapidity. In Sec. I VII we 
discuss the results with the running coupling, both for the case without and with impact parameter dependence, and 
for two different prescriptions of the running coupling. Finally, in Sec. I VII] we state our conclusions. 



II. BK KERNEL IN LO AND BEYOND 

In the leading logarithmic approximation in Inl/a:, the nonlinear Balitsky-Kovchegov [H, [l^, [l8l - [2^ evolution 
equation derived in dipole picture [l7j has the following form 

= "^y ■^(xo-x2)2(xi-x,)2 + - - ' 

where as — agNc/n is the strong coupling constant. Here, A'^xoxi = Af(xo,Xi,y) is the dipole-nucleus scattering 
amplitude, and Xq, xi are two-dimensional vectors of the transverse position of the dipole ends. Alternatively, one can 
introduce the vector denoting the dipole size rgi = Xg — Xi, and the impact parameter boi = (xq + Xi)/2. Thus in 
general, the amplitude depends on the four degrees of freedom in transverse space and rapidity, Y ~ ln(l/a;), playing 
the role of the evolution parameter. The transverse part of the LL kernel 

dz (Px2 Xqi 

z 27r a;o2'^i2 

is conformally (Mobius) invariant in 2-dimensions. Here, we introduced a more compact notation denoting x^ = 
Xi — Xj, Xij — |xy I and z is the longitudinal momentum fraction so that rapidity is y — In I /z. 

To obtain the solution of this equation, one has to specify an initial condition aXY = Yq: A^xoxi = A^'-^^xo, xi; Y — 
Yo). 

The amplitude N^^xi in ([Ij is given by the following correlator 

^x„x, = -^Tr(l - C/t(xo)C/(xi)) , (2) 
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where the trace is done in the color space, and the eikonal factor U is defined as the path ordered exponential with 
the SU{N) gauge fields (in the gauge A~ = 0) 

C/(x) = Pcxpl^iJ da;-T''A+(a;-,x)| . (3) 

The averaging (...) in ([2]) is performed over an ensemble of classical gauge fields. 

The next-to- leading logarithmic corrections to the BK equation have been computed in [13, H^. It would be 
instructive to perform the analysis of the BK equation with full NLL approximation. Due to the complicated form 
of this kernel we are not trying to solve it here. What we consider instead, is the equation with two important 
modifications. The first one is the running coupling correction, which we analyze in Sec. IVII of this paper. The second 
one is the subleading correction coming from the kinematical corrections, included in the form of the modified kernel 
proposed in (43| . 

The soft gluon approximation performed in [l7j was crucial in order to factorize the the branching kernel for 
the gluon emissions. In [4^ a modified version of the kernel was proposed, which includes part of the subleading 
corrections coming from the kinematics. By taking into account an improvement in kernel in which the value of 
energy denominator is obtained from the invariant mass of the produced gluon pair. This leads to a new color dipole 
kernel which has the form 



dz d X2 
— z — 7^ — 



X02 



Xl2 
Xoi 



2Ki 



X02 

xoi 



Xl2 
Xoi 



X02 ■ Xi2 



(4) 



where Ki is the modified Bessel function. The modified kernel becomes equivalent to the leading logarithmic kernel 



for zxq2 ^ x"^ 



01 



and 



t'21 



^ Xqi but differs significantly otherwise. To be precise, the production of larger dipole sizes 



is exponentially suppressed above the cutoff size which depends on the longitudinal momentum fraction of the soft 
gluon. The modified kernel introduces corrections at all orders beyond the leading logarithmic approximation. The 
suppression of the large dipole emission implies also a suppression of the diffusion in the impact parameter. An alert 



reader might notice that the above modified kernel is not conformally invariant. This is due to the fact that in [43 1 
only part of the kinematical corrections were taken into account resulting in kernel (j4|). In general there should be 
also a cutoff on small dipoles, similarly to what was done in [50| . which would yield conformally invariant kernel. We 
postpone the derivation and numerical analysis of such kernel to a later work. 

In what follows we will analyze both forms of the equations numerically, taking into account also running coupling 
corrections. We shall demonstrate that the sensitivity to the subleading corrections strongly depends on the type of 
approximation one is working on, i.e. whether the solutions are impact parameter dependent or not. 



III. NUMERICAL METHODS 



The BK equation was solved numerically by discretizing the scattering amplitude in terms of variables 
(logj^o logio ^' '^o^ where 9 is the angle between impact parameter b and dipole size r. There is an additional 
angle, which refiects the overall orientation, we assumed that the solution is rotationally invariant. The amplitude 
N{r, b, cos9) was placed on a grid with dimensions 200^ x 200;, x 20e. It was found that the number of points in the 
grid was less important than the absolute limits on the grid and the number of integration regions. The limits of the 
grid were 10~^ — >■ 10^ in r and b. We note that, the grid limits of r and b cannot not be set independently due to 
the correlations in dipole size and impact parameter. The initial condition A^*^"^ at rapidity Y = was chosen as in 
Ref. jU] to be of the Glauber- Mueller [HI, IHi] form 



^r(0) = l_g-10r^e-V ^ (5) 

This is the form of the scattering amplitude at the initial rapidity. We note that in this analysis the constants are 
arbitrary and have not been fit to any experimental data. The initial condition ([5]) has the property that it reaches 
unity when the dipole size is large and it goes to zero for large values of impact parameter by the steeply falling profile 
in b. The amplitude is evolved in small steps Ay = 0.2 as it was found that the smaller spacings had no significant 
effect on the accuracy of the evolution of the amplitude. Below, we denote integration over X2 simply by J^^ and this 

is understood to absorb the kernel as well as a measure d^X2 for ease of notation in this section. The basis of the 
method for the solution is the same as utilized in Ref. [4l| and we briefly outline it below. After one step of evolution 
from Yq to Y = Yo + AY the equation ([T|) becomes 
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iV,oxi(>o + Ay) = Ar^o^i(^o) + / / [iV,„,, + N,,,, - N,,,^ - N,,,,N,,.,,] . (6) 

The lowest approximation for this equation is 



X2 JYa 



^(1) = ^(0) ^ / [A^(") + iV(°) ~ iV(°) - A^f") iV(°) 1 (7) 

a:o2;i xqx\ ' / L-'*a:oa:2 x\X2 xnx^ xqX2 x\X2\ V'/ 



a:o^l ^03^1 / L a:oa:2 a;ia:2 ^o^^l Xi^X2 x\X2\ 

■J X2 

The superscript 1 denotes the first approximation while the superscript denotes the initial condition. The accuracy 

of the equation is measured by the relative difference of the approximation defined by | (N^c^Xj — Nx'^l- ) /Ni^l- \ . If this 

is less than some required accuracy then the value of Nx]lj is taken as the scattering amplitude at F = Iq + Ay. 
In our case the relative accuracy was set to be equal .02. If the accuracy condition is not satisfied, then we can 

find a second approximation by utilizing linear interpolation N^ll. = — ^ — ""'"af "'^"^'^ -^^i^j ■ using this linear 
interpolation in the right hand side of [6] the expression for the next iteration can be attained 



^(2) ^ ^(0) I ^ /" mi^) I 7V(0) _ ^(0) _ ^(0) ^(0) 

XqX\ ^^XqX-i ' rj / V X(iX2 ' X\X2 XqX-i XqX2 XlX: 

'X2 



2 

Ay 

~2~ 



[iV(i) +iV(i) -7V(i) -7V(i) 7V(i) 

L^'xoX2 ' XiX2 ^^XqXi XoX2 XiX; 



+ ^ I Wgl~Ni%jN(\l^Ni%J. (8) 



X2 



If the accuracy condition mentioned earlier is still not satisfied this can be repeated by utilizing Af(") = 

HE. — ^y~^ — ^ + N^^K The number of iterations required to reach this level of accuracy varies on how close one 
is to the initial condition, and the exact form of the equation (kernel used, running coupling form, etc). Near the 
initial conditions it could take upwards of ten iterations but far away from the initial condition two or three iterations 
are usually sufficient to obtain the accuracy assumed. With this in mind an algorithm was written so that the iteration 
number was not fixed but only continued until the desired accuracy was reached. Once the amplitude has been found 
for y = yo + Ay that takes the place of A^(°) in the numerics and the procedure is repeated to find the amplitude at 
y = yo + 2Ay. in this way the amplitude is evolved numerically to higher rapidities. 

As was mentioned earlier the number of points in the grid was not as important as the number of integration 
regions the program executes to get a good solution. For our limits it was found, that breaking the range into at 
least 20 integrations was needed. This meant there would be an incredible number of function evaluations and in 
order to reach any useful rapidity on a standard dualcore machine would take upwards of a month. The program was 
parallelized using the MPI (Message Passing Interface) library so it could be used with computer clusters to run the 
program. Each run took approximately three days to run on 32 linked 3 GHz processors and generated a total of 5 
gigabytes of data. 



IV. RESULTS WITHOUT IMPACT PARAMETER DEPENDENCE 



In this section we briefiy present the results of the numerical solution without impact parameter dependence. This 
was of course done in many previous works [ssl - lssl . [ssj . but we repeat this analysis below for several reasons. First of 
all, we wanted to have benchmark for the comparison of the solutions. Secondly, we wanted to perform the analysis 
of the modified BK equation with Bessel function kernel, a computation not performed earlier in the literature. 

The BK equation can be evaluated without impact parameter or angular dependence and this greatly reduces the 
time needed to evolve the scattering amplitude. The initial condition used in this case is taken to be 



7V(") = 1 - e-"' , (9) 

without an impact parameter profile. This form goes to unity as the dipole size gets large and has a narrow transition 
region from to 1. The absolute limits on the dipole size here are lO^^'^ — > 10^'^ and the number of points in the grid 
is 5000. These results are for fixed coupling = 0.1, the case with running coupling is considered in section Sec. IVII 
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FIG. 1: Left: Comparison between the LO kernel (solid lines) and the kernel with Bessel functions, ([4]) (dashed lines) at 
constant rapidities. The solid line at the far right is the initial condition for both cases and each line to the left represents an 
evolution in ten units of rapidity to a maximum of Y = 40. 

Right: plot of the logarithm of the saturation scale as a function of rapidity. The solid line represents the evaluation with LO 
kernel and the dashed line represents the solution with Bessel kernel. The slope of the curve gives As and the three regimes of 
evolution can be seen. Calculation done for the fixed coupling case and without the impact parameter. 



The solid lines shown in left graph in Fig. [T] correspond to a set of scattering amplitudes found from the LO BK 
equation at constant rapidity. The result has the well known form of a traveling wave with rapidity playing the role 
of the time in the evolution. As rapidity increases the front moves towards the small values of dipoles. 

Next, we performed the evolution of the modified equation with the Bessel function kernel and the solutions are 
illustrated by the dashed lines in Fig. [T] The normalization of the amplitude in the dilute regime is much smaller 
than the LO case. This is to be expected as the Ki functions fall off faster than the power-like LO kernel. The speed 
of this evolution can be quantified by the evaluation of the saturation scale. The equation to define saturation scale 
is 



(iV(r = l/Q„y)) = AC, (10) 

where k is a constant between and 1. It determines the point at which nonlinearities begin to become important 
in the evolution equation. In all the analysis in this paper n = 0.5 was chosen. Based on results of [54|| we expect the 
saturation scale to have the following form 



Ql{Y) ^ Ql exp 



Xilc) 

l-7c 



Y 



2(1 -7c) 



InF 



(11) 



where x{l) is the BFKL kernel eigenvalue, and 7^ — 0.62. This gives 



X(7c 



and 



2.39 for the LO 



l-7c 2(l-7c 

kernel. The second term in the exponent is a correction of less than 10 percent compared to the first term for the 
rapidities which we are considering so we take the saturation scale to be parameterized in this paper by 



where As is extracted from the numerical solution and Qo is a normalization term. The extracted value for the 
saturation exponent is As — 4.4 for both the LO and Bessel function kernels. The exponent was extracted both for 
7xs = 0.1 and cfs = 0.2 and it was found to be the same. The values of As are evaluated in the region of rapidities 
where the exponent is constant, i.e. nearly asymptotic. We observed that for the lower rapidities, near the initial 
condition the exponent slowly increases and at very large rapidities the exponent begins to decrease. The latter effect 
is due to the absolute limits of the dipole size (or the box in which the equation is being solved). This strongly affects 
the solution, and one of the effects is to slow down the evolution. This effect can be clearly seen in Fig. [1] where 
the three regimes are distinct. The first regime is where there are a lot of preasymptotic effects and there is strong 
dependence on the initial condition and the speed of evolution increases (approximately between and 34 units of 
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rapidity), the regime where evolution is a constant (approximately between 34 and 47 units of rapidity) and finally 
the regime where box size limits growth (rapidity 47+). These values clearly shift with box size and they also change 
depending on the value of and to a lesser extent which kernel is used. 

The extracted value of the exponent of the saturation scale is consistent with the theoretical predictions of 'b^. 
What is however surprising at first, is the fact that the asymptotic value of is nearly the same for the modified 
version of the equation. This is also evident by investigating right plot in Fig. [1] where the two saturation scales grow 
with rapidity in almost exactly the same way. This effect can be easily understood by inspecting the form of the 
modified kernel 



Qoi 



Kf{QoiXo2) + Kf{QoiXi2) - 2Ki{QoiXo2)Ki{Qoixi2) i^"^,,^^^, 

F02IF12I 



(13) 



Here Qoi = y/z/xoi and z = e~^. When Qoi is small and the arguments of the Bessel functions are small this 
expression reduces to the LO kernel. At large rapidities where the exponent Ag is evaluated, z is small and thus Qoi 
is also small. The difference between the modified kernel and the LO kernel comes from the regime where the dipole 
sizes a:i2,a;o2 are large, that is when 

1 

a;i2,a;o2 > 7;— • 

WOl 

In this large dipole regime (and for large rapidities), the scattering amplitude is close to unity and the right hand 
side of the BK equation goes to zero so there is no contribution from that regime. This means that the effect of 
the infrared modification of the kernel of the type is negligible when the nonlinear evolution is considered. The 
modified kernel (fT3)) was shown to reproduce the double-logarithmic terms in the exact NLL calculation. These type 
of terms are known to come from the scale choice in the amplitude, (55l - [57| . 

Two comments are in order here. The first observation is that, as we shall see later, the presence of the impact 
parameter changes the dynamics significantly, and in this case the dependence on the form of the kernel is more 
pronounced. In particular, in this case we do find the differences in the evolution speed due to the cutoffs in the 
infrared. The second comment is that, in order to address the question of the scale choices one needs to take into 
account the corrections beyond the BK equation. BK equation is highly asymmetric with respect to the target and 
projectile and therefore cannot address the problems of all the scale choices. The corrections would necessarily have 
to include the Pomeron loops in order to make the evolution symmetric with respect to the target and projectile. 
Such formulation has been implemented in a Monte Carlo approach [isl . Isol . Issl . [s^l with Pomeron loops effectively 
taken via so-called dipole swing and with cutoffs both in the infrared and ultraviolet regions. 



V. RESULTS WITH IMPACT PARAMETER DEPENDENCE 



A. Dependence on the dipole size 



In this section we discuss the results which include the impact parameter dependence as well as the angular 
dependence between the dipole sizes and the impact parameter orientation. Let us first investigate the dependence 
of the scattering amplitude on the dipole size with fixed impact parameter value. Fig. [2] shows that case of the LO 
simulation with ois — 0.1. The scattering amplitude at small dipole sizes evolves in a similar manner as in the case 
without the impact parameter. 

At large dipole sizes the situation is drastically changed with respect to the translationally invariant case. Here, we 
observe that the amplitude drops down from the initial distribution and forms a second evolution front. This drop 
is very rapid compared to the expansion evolution of the small dipole size regime or the large dipole size regime at 
higher rapidities. The evolution of the front at large dipole sizes can be best seen in Fig. |2(b)| where the steps in 
rapidity are greater. As discussed in [4l| there is a clear physical reason for this effect. For a large dipole, its end 
points are in the region where there is no gauge field. In this situation the gauge field correlator 



Ar(x,y) = i-Tr(l-f/t(x)f/(y)) , (14) 

vanishes because C/(x) ~ U{y) ~ 1. In this case the dipole is larger than any of the other scales in the problem, 
including the impact parameter of the collision. In other words this situation corresponds to the setup of very large 
dipole scattering on a localized target, and in this case the dipole misses the target, resulting in the vanishing scattering 
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Impact parameter: 1 .000 | cos{e): 0.0 | AY: 1.0 | max Y: 10.0 



Impact parameter: 1 .000 | cos(9): 0.0 | AY: 1 0.0 | max Y: 50.0 





10 10 10' 

Dipole Size 



lO^IO'^lo'^IO'^IO^'lO'^IO'^IO 10° 10' 10^ 10^ lo" 10^ lo'' 10^ lo" 
Dipole Size 



FIG. 2: Graphs of the scattering amplitude N{x, b, 6, Y) at various constant rapidities for the calculation with impact parameter 
dependence. Left: the individual lines correspond to the rapidity intervals AY — 1.0 up to rapidity of K = 10. Right: the 
individual lines correspond to the rapidity intervals of AY = 10 and the maximal rapidity is K = 50. The dotted line represents 
the initial condition which is the same in both graphs. Impact parameter and angle are fixed for both graphs at b = 1.0 and 
cos(6') = 0. 



Impact parameter: 0.110 | Gos(e): 0.0 1 AY: 5.0 | max Y: 50.0 



Impact parameter: 1 00.000 | cos(e): 0.0 | AY: 5.0 | max Y: 50.0 




Dipole Size 
(a) b = 0.11 




'iV 10"^ 10"* 10"^ 10'^ 10'' 10° 10' 10^ 10' 10* 10^ lO*" 
Dipole Size 



(b) b = 100 



FIG. 3: Graphs of the scattering amplitude at different constant rapidities for two different impact parameter values. The 
dotted line represents the initial condition which is the same in both graphs, however the initial condition is near zero on the 
second graph with cos{6) = 0. Each line past the initial condition represents five rapidity units to a maximum of 50. 



amplitude. This effect is not present in the case where impact parameter dependence is neglected simply because this 
approximate case corresponds to the infinitely large target. 

It is interesting to investigate the dipole size dependence of the amplitude at large values of the impact parameter. 
There the initial condition, sets the amplitude to zero. The evolution quickly changes this value, and a unique feature 
of the solution develops. Namely, a peak is formed with the center at the dipole size value which is exactly twice the 
impact parameter value. The peak grows until saturation is reached and the evolution of the fronts proceeds to the 
infrared and ultraviolet regions. The discussion about the origin of peaks is given in subsection IV Dl 

We have repeated this analysis for the case of the modified kernel. As can be seen in Fig. 4(a) the evolution in the 
small dipole regime is now quite different from the LO kernel. The evolution is significantly slower in this regime for 
the Bessel function kernel. Interestingly, the evolution in the large dipole region is slowed too but not as much as for 
small dipoles. This is best illustrated in Fig. |4(b)| where the evolution has been performed to large rapidities. The 
origins for this behavior are discussed in detail in the section on the saturation scale and evolution speed. 
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Impact parameter: 1 .000 | cos(e): 0.0 | AY: 2.0 | max Y: 10.0 Impact parameter: 1.000 | cos(9): 0.0 j AY: 10.0 j max Y: 50.0 

— — — — — I I I I I I I I I I I 
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FIG. 4: Graphs of the scattering amphtude as a function of the dipole size at various constant rapidities for fixed impact 
parameter b — 1.0 and angle cos{9) = 0. Solid lines are for the LO kernel and the dashed lines correspond to the Bessel kernel. 
The initial distribution is equivalent for both kernels and is represented by the dotted-dashed line. On the left graph each line 
represents a change in two units of rapidity to a maximum of ten and on the right graph each line represents a change in ten 
units of rapidity to a maximum of fifty. 

Dipole Size: 1.000 | cos(9): 0.0 | AY: 10.0 | max Y: 50.0 Dipole Size: 0.110 | cos(e): 0.0 | AY: 5.0 | max Y: 30.0 
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FIG. 5; Graph of scattering amplitude as a function of impact parameter for fixed dipole size r — 1.0. The solution with the 
case of the LO kernel is plotted as a solid line and with the modified kernel (|13p as a dotted line. The dotted-dashed line on 
the left is the initial condition. Each line thereafter represents an increase in rapidity of ten units to a maximum of fifty. Right 
plot: the same but for the dipole size r — 0.11 and in logarithmic scale for the amplitude. 

B. Impact parameter profile of the scattering amplitude 

Dependence of the dipole amplitude on the impact parameter is illustrated in Fig. [5l The leftmost dashed-dotted 
line is the initial condition Eq. [S] which has a very steep profile in impact parameter. The evokition of the scattering 
amplitude towards large values of impact parameter follows the diffusion of large dipoles. The speed of this evolution 
can be extracted numerically and is determined by the expansion of the black disc radius. We will discuss this quantity 
in detail in the next section. 

Evolution in impact parameter shows a marked change in profile from the steeply falling exponential in the initial 
condition. This is better illustrated in right plot in Fig. [5]where we replot the impact parameter using the logarithmic 
scale in scattering amplitude. The profile changes from the exponential to a power tail at small scattering amplitudes. 
This can be seen as an 'ankle' in the curves of constant rapidity. The origin of this power-like tail was discussed 
in detail in Ref. [4lj. These power tails are also present in the modified kernel. In the latter case however there 
is a slower evolution of the profile towards the large values of impact parameters. There also exists a nontrivial 
angular dependence which is most prominent in the cases of large dipole size or impact parameter but for very specific 



FIG. 6: Two configurations, 'aligned' and 'perpendicular', of tlie dipole orientation witli respect to the impact parameter. 



Impact parameter; 100.000 | cos(e): 1.0 , -1.0 | AY: 10.0 | max Y: 50.0 Impact parameter; 100.000 | cos(9); 0.0 | AY: 10.0 | max Y: 50.0 
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(a) cos(6») = 1.0, -1.0 (b) cos(e) = 0.0 

FIG. 7: Graphs of scattering amplitude versus dipole size for fixed impact parameter b — 100.0 and various rapidities and 
angles. The initial condition is the same in all graphs and it is near zero, each curve represents an increase in ten units of 
rapidity to a maximum of fifty. The LO kernel (solid lines) and the Bessel kernel (dotted lines) are plotted on the same graph. 

configurations. In the case when the dipole size is much smaller or much larger than the impact parameter the solution 
does not depend much on the spatial orientation of the dipoles. On the other hand, for the case when the dipole size 
is twice as large as the impact parameter there exists strong angular dependence. These effects are best illustrated in 
Figs. [7Ua,b) and[5Ja,b). The scattering amplitude as a function of both dipole size and impact parameter for different 
choices of the angle 9 is shown, where 9 is defined as the angle between the dipole and the impact parameter, as 
illustrated in Fig. [SI The amplitude has a peak when at r = 26 for both orientations as is shown in both plots in Figs. [T] 
Note however, that the peak is distinctively sharper for the 'aligned' dipole configurations, when cos9 = 1.0 or — 1.0, 
than for the 'perpendicular' configuration. This is also illustrated in Fig. |9] where the dependence on the angle is 
shown. For values that are near the r — 2b point there are enhancements at cos 6* = 1.0, —1.0 and this is present in 
both kernels. These effects can be seen in both plots in terms of dipole size and impact parameter. It is interesting to 
note that the peak is present in the case of scattering amplitude versus dipole size the peak even when cos 9 = 0.0. On 
the other hand such structure is absent for this configuration in the impact parameter profile with fixed dipole size. 
It is also evident in Fig. [5] that the amplitude is fiat in impact parameter when the dipole size is much larger than b. 
We will demonstrate that all these effects can be easily understood from conformal representation of the amplitude 

C. Saturation Scales 



The saturation scale in the impact parameter dependent scenario is again defined by the following equation 

{N{r^l/Qs,b,9,Y)) = K, (15) 
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FIG. 8: Graphs of scattering amplitude versus dipole impact parameter for constant dipole size r — 100.0 and various rapidities. 
The initial condition is the same in all graphs and it is the steeply falling dotted-dashed curve, which is the same for both the 
evolution with LO kernel (solid lines) and the Bessel kernel (dotted lines). Each curve represents an increase in ten units of 
rapidity to a maximum of fifty. 



Dipole Size: 1.000 I impact parameter: 1.000 | AY: 10.0 | max Y: 50.0 Dipoie Size: 100.000 | impact parameter: 47.863 | AY: 10.0 | max Y: 50.0 




-1.0 -0.5 0.0 0.5 1.0 Cos(e) 

Cos(e) 

(b)r = 100.0 b = 47.8 

(a)r = 1.0 f) = 1.0 

FIG. 9: Graphs of scattering amplitude versus angle. Both have the LO kernel (solid line) and the Bessel kernel (dashed 
line) graphed on them. The first graph shows no angular dependence and the second shows a marked increase at values of 
cos{9) — 1.0, —1.0 which accounts for the difference in peaks in Figs. (|8I7|) 



where k is a constant. In all the following analysis we have set k = 0.5. It is important to note that, in this case 
the form of the amplitude admits two solutions to the above equation. As is evident from Fig. [5] one solution for the 
saturation scale is for a larger dipole size and one for a smaller dipole size. The saturation scale Qs always refers 
to the solution where the dipole size is smaller. We have found that the slope in rapidity of the saturation scale 
Qs increases for low values of rapidities, then reaches an approximately constant value and for ultrahigh rapidities it 
starts to decrease. The first effect is caused by the preasymptotic contributions, the latter effect is caused by the finite 
size of the grid. We have found that the effects of the grid can be neglected below the rapidities of order ^ 60. The 
saturation scale as a function of the rapidity is shown in left plot in Fig. (TU] The solid line shows the calculation in 
the case of the LO kernel and the dashed line is for the Bessel kernel. It is clear that, the dependence on the rapidity 
is exponential as expected for the computation with fixed value of the coupling. The numerical value of the exponent 
governing the rapidity dependence of the saturation scale, extracted in the LO kernel case is = 4.4, see pB)) . The 
value of the exponent extracted for the evolution with Bessel kernel was found to be As = 3.6. Clearly, the sublcading 
effects of the modified kernel cannot be neglected here, which has to be contrasted with the case without the impact 
parameter. 
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FIG. 10: Left: plot of the saturation scale as a function of rapidity. The solid line corresponds to the LO simulation and the 
dashed line is the saturation scale of the Bessel kernel. The impact parameter is fixed at 6 = 1.0. Right: the dependence of 
the saturation scale on impact parameter for two different values of rapidity. Solid line: LO kernel, dashed line Bessel kernel. 
Strong coupling is fixed at cfa =0.1. 



In the case when impact parameter is much larger than the inverse of the saturation scale the exponent in rapidity 
is independent of the impact parameter value. This means that the saturation scale has a factorized form 

Q^Y, h) = Ql exp(a,A,r)5(6) . (16) 

This is demonstrated in Fig. [TU] where the small dipole saturation scale is shown as a function of the impact 
parameter for two different values of rapidity. The power tail 1/6'* is clearly prominent. There is significant difference 
between the saturation scale from the Bessel function kernel and the LO kernel. 

The second solution of the equation ([TSl) . which shall be called Qsl{Y, b), gives saturation scale at large dipole size. 
The same analysis that was performed on can be performed on Qsl with the parameterization of this saturation 
scale taken to be 



QL = QoLe-'=^^=^, (17) 

where once again Qol is a normalization term and the minus sign in the exponent is because the evolution is now 
moving towards larger dipole sizes. 

The extracted value for the LO kernel is Xsl — 6.0 and for the Bessel function kernel Xsl = 5.6. The difference 
between these two exponents is now about 7% and it can be seen on Fig. [TT] that these curves are much closer than 

in Fig. nni 

The reason for these effects can be understood by again inspecting the form of the Bessel kernel 



Qoi 



Kf{QoiXo2) + Kf{Qoixi2) - 2Ki{Qoixo2)Ki{QoiXi2) 



X02 ■ Xi2 

|a;o2||a:i2| 



(18) 



When 3:01 is small then Qoi is large and the cutoff is on dipoles such that 0:02 > 1/Qoi is very large. This means 
that unless Xq2 is very large (which would also correspond to large X12) then the modified kernel is very close the LO 
kernel. It is this limitation of the phase space of the modified kernel that causes the dipole saturation scale for large 
dipoles to have a rather modest difference between the Bessel and LO case. At this point it is worth to note that the 
Bessel kernel does not exhaust all the kinematical effects. To be more precise, we would expect that the kinematical 
cuts are resulting in the kernel which is also conformally invariant. This will give cuts on both small and large dipole 
sizes and further reduce the evolution speed. 

One can also define from (jlSp a scale which corresponds to the extension in impact parameter space. This scale is 
the radius of the black disc. This can be done by solving this equation for b rather than for the dipole size. 

We therefore define the black disc radius in impact parameter Bs by solving the equation 
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FIG. 11: Left: plot of the saturation scale for large dipole size as a function of rapidity. The solid line represents the result 
using LO kernel and the dashed line represents the result using the Bessel kernel. The impact parameter is fixed at fe = 1.0. 
Right: plot of the black disc radius defined by Eq. [l9]for the LO kernel (solid line) and the modified kernel (dashed line). The 
dipole size is fixed at r = 1.0 
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TABLE I: Summary of extracted saturation exponents for solutions with impact parameter. (1) means Glauber-Mueller initial 
conditions, (2) means Glauber -Mueller with the exponential cutoff on the large dipole sizes. 



{N{r,b^B,,e))^K, (19) 

with respect to b, and where once again k = 0.5 is chosen. We assume the exponential form for the behavior of the 
impact parameter radius as a function of rapidity 

^ 52^6^-^^=^ . (20) 

Here B^q is a normalization term and Xbs is extracted from the numerical solution. We have found for the LO 
kernel Xbs = 2.6 and for the modified kernel A^.; = 2.2. This is approximately half of the As in the case of the LO 
kernel. As argued before [111, [60| this is due to the fact that the amplitude depends on one variable, which in the case 
of the configuration 6 3> r is proportional to r/b'^. This immediately means that the impact parameter dependence 
in rapidity is twice slower than the one of the dipole size. The computation for the Bessel function kernel shows that 
it does not hold as closely in this case because in such case we do not have an exact conformal symmetry. These 
properties will be explained in more detail in the next section. 



The simulations with different value of the fixed coupling were also performed. The results are summarized in table 
m It can be seen that the exponents for the LO kernel do not depend on the value of the coupling, which is the 
expected behavior as this is a fixed order calculation. On the other hand the exponents extracted for the modified 
kernel significantly differ exhibiting the nonlinearity in the coupling constant. The exponents are further reduced with 
respect to the LO values which is due to the resummation of the subleading terms in \nl/x. 

The dependence on the initial conditions was also tested. As an alternative, we have taken the second initial 
condition to be modified by the cutoff in the large dipole size 

Afg = Af^°' exp(-r^) , (21) 
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(a)Impact parameter dependent amplitude for the fixed coupUng 
as = 0.1 case as a function of the dipole size and for different 
values of the impact parameter. 



(b)Saturation region in (Y, Inr) plane for different impact 
parameters. 



FIG. 12: 



•where /j, = 1/5. The exponents are also sho-wn in table U (LO (2)). We observe a modest variation of the exponents 
with the change on the initial condition. The most significant change is for the large dipole saturation scale. This is 
to be expected as in this region the two initial conditions differ significantly. 

In Fig. ll2f a) the amplitude as a function of the dipole size and various impact parameters is shown. It is interesting 
that for large dipole sizes the amplitude has the same front for all the impact parameters. This is related to the 
properties of the solutions stemming from the conformal symmetry, see Sec. lVDl In Fig.[T2lb) we show the saturation 
region for the solution with impact parameter. Unlike in the local case (without the impact parameter) here the 
saturation region has a 'V' shape in (y,lnr) space, which is moving towards higher rapidities and larger dipole sizes 
as the impact parameter increases. Different shaded areas correspond to three different impact parameters. Again, the 
common front for the different values of b is clear. The distortion at lower rapidities and for small impact parameter 
stems from the initial conditions. 



D. Conformal representation and properties of the amplitude. 



Most of the features observed in the numerical solutions can be explained by using the conformal representation of 
the solution for the scattering amplitude. In general the representation can be shown to be of the form Q 

rc+QO 1 

F(ro,r,b;y) = / — cxp(c^r) F(ro, r, b; w) , (22) 



with 

-f oo 



F(ro,r,b;a;) = ^ / dv I £^ iv"^ 



4 J p + (n^)2] ^ + ("±1)2] 

X E--* - ^' -y - ^) ^"''^ (I + b - w, -| + b - w) , (23) 

where Tq, r are the transverse sizes of two scattering objects (for example onia), b is their relative impact parameter. 
The conformal eigenfunctions are defined as 

i?"-(pio,P2o) = (-1)" , (24) 

\P10P20J \PiqP20/ 

where complex notation for the two dimensional vectors {p^, Py) has been used 

p = Px + ipy, P* = Px- ipy , 
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and where the conformal weights are 



1 — n - 1 — 71 
h = — \- iv, h — — h w . 



Function J-{i'^ n, a^; uj) contains the details of the dynamics. For the case of evolution with linear BFKL, the form of 
it is well known 



-BFKL^ - 1 



u} - x{n, v) 



with 



being the LO BFKL kernel eigenvalue. In the case of the nonlinear equation the exact form of the function 
J-{v,n^as',oj) is unknown. The origins of the peaks in the amplitude can be understood by analyzing the trans- 
verse structure encoded in functions E^-'^ . We fix rg and investigate the dependence on r from the transverse integral. 
We switch from the vector notation to the complex notation for the arguments of the E functions. Using the explicit 
expression ((24)l we obtain 



where we switched to the complex notation for the arguments. The biggest contribution comes from the region of 
w ~ 0. For our purposes it is also enough to take n = 0. In this region the integrand has the approximate form 



Using b = |5|e'^', r = \r\e''^'- and AO = Or - Ob we have that 



1/2+ii' 

rr^ ^1 ■ (27) 

It is immediately clear that there will be angular dependence for the & — § case with the configurations of 'aligned' 
dipoles A9 = 0,tt giving largest contributions. In the case of the 'perpendicular' orientation of dipoles with respect 
to the impact parameter 9 — 7r/2, — 7r/2, the expression reduces to 

(^)"''"- 

This structure is responsible for the presence of the peak in the amplitude in the case when the b is fixed and r varied, 
and the absence of the peak in the case when r is fixed and b varied, for A6 = 7r/2, —t:/2. This corresponds to the 
situations in right hand plots in Figs. [7] and H] correspondingly. 

An expression for the saturation scale dependent on the impact parameter can be derived using the method in 
[ssj . To this aim one needs to take the Mellin representation for the solution to the linear equation and apply the 
absorptive boundary. The integral over the transverse variable can be performed using the representation [j] 

d^p^E'^'^ipio, P2o)-B"''^*(pi.o, P2'o) = cix''x*'^F{h, h, 2h; x) F(h, h, 2h; x*) 

-I- C2X^~''x*^~'^F{l - h,l- h,2-2h; x)F{l - h,l - h,2 ~ 2h; x*) , (29) 
where 

P12PV2' .„„N 

X = , (30) 

Pll'P22' 
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is the anharmonic ratio and F are the hypergeometric functions. To obtain the saturation scale we take n = 0, 
and expand around a; ~ 0. This simphfies the above expressions as in this hmit F(h, h, 2h; x) ^ 1 and the whole 
dependence on x comes through factors x^x*^. In the case when the impact parameter b is much larger than the 
dipole sizes, 6 » r,rQ, one has 

N = ^. 

Putting everything together, the scattering the amplitude in the linear evolution case reduces to [6l|, [g^] 



F(ro,r,b;r) = -%^ e^p ( a^xWY + 2^l.\n ^ ] . (31) 



Here, we have taken into account the contribution from only zero conformal spin. Using A = 1/2 + w the above 
expression can be recast into 

F(ro, r, b; Y) = -a^ ^ 2(A - 1/2) exp U,x(A)r - (1 - A) In j j , (32) 

where the prefactor in pip has been expanded around v = Taking the saddle point condition and the condition 
that the exponent vanishes at the saddle point which is the requirement on the saturation boundary one arrives at 
two conditions for this line (noted by a subscript 0). 

asY^{\^)+M\?l{n,r)f = 0, (33) 
a,rx(Ao)-(l-Ao)ln(6V(ror))2 = 0. (34) 

These equations can be solved to yield the saturation scale but it was found that one can include further corrections 
[11]. We can obtain these corrections by using the solution to the saddle point equation to find Ac, and use this to 
evaluate the prefactor. The resulting modified equations are then 

a3rx'(Ac)+ln(6V(^oO)' = 0, (35) 
a,rx(Ac) - (1 - Ae) \x,{}?IMf - ^ ln[a,yx"(Ac)] • (36) 

By keeping one of the dipole sizes fixed, say tq, we can solve for r to get the saturation line 

Q2 (ro, h- Y) = g 3 . (37) 

For large cisY the Ac approaches Aq value, with Aq = 0.37. The saturation scale has 1/6'* dependence which comes 
automatically from conformal symmetry. One can solve the above equation for h and keep tq, r fixed which yields 



e 2(1-Ac) 

S2(ro,r;y) = r^r . (38) 

[a^y'x"(Ac)]*<i^^=) 

This is the rate of the expansion of the radius in impact parameter space. Note that, the speed of the expansion is 
governed by the exponent which is half that of the saturation scale and the dependence on the dipole size is linear. 
This is also found in the numerical solution. For the large dipole sizes r 3> tq, 6 the anharmonic ratio reduces to 

N = ^^. 

r 

Following the same scheme one obtains for the saturation scale 

Qh (ro ,b;Y)= ^ , . (39) 

r2[a.rx"(Ac)]^^™ 
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Modified Kernel = 0.1 


2.4 
2.0 


LO Kernel cts = 0.2 
Modified Kernel Us = 0.2 


2.6 
1.6 



TABLE IL Extracted exponents governing the behavior of the black disc cross section. 



The saturation scale for large dipole sizes is independent of the impact parameter b. This is also found in the solution, 
as is clear in Fig. 1121 Therefore the 'V shape of the saturation region is a consequence of the conformal symmetry of 
LO kernel. From the above considerations one can see that the rapidity behavior of both saturation scales is identical 
for large and small dipoles. We found that the two exponents differ somewhat, see TableU) Most probably this is due 
to the initial conditions which are asymmetric in both small and large dipole sizes. However, a more detailed analysis 
is needed to confirm this effect. 

In general we see that, both saturation scales, are in fact originating from one saturation scale due to the fact that 
the solution is expressed in terms of the anharmonic ratio. 



E. Dipole cross section and blaclc disc radius 



By integrating the amplitude over the impact parameter the dipole cross section is obtained as a function of the 
dipole size and rapidity. Despite the fact that the amplitude is bound and never exceeds unity, the dipole cross section 
can still increase very fast due to the fact that the amplitude has power tails in impact parameter, see discussion 
in [6^, 64] . We thus expect the power like growth of the dipole cross section with the energy, or exponential with 
rapidity. 

The dipole cross section is defined as an integral over b of the amplitude 



cr(r,y) ^ 2 J (fhN{r,h,Y) . (40) 

In what follows, we will investigate the part of the dipole cross section which is coming from the black disc regime. 
To be precise, we integrate the amplitude over the values which are close to unity. This is once again performed by 
constraining the amplitude through the equation p9)) . 

The black disc part of the cross section is therefore defined as 



TBBir^Y) = 2 J <fhN{r,h,Y)e[Nir,h,Y) - K,] « 2nRlj,{r,Y) , (41) 

This black disc cross section is plotted in Fig. [T3] and it can be seen that the slope of the black disc cross section 
in rapidity reaches a constant value at large rapidities. One can parametrize Rbd as 



Ri^{r,Y) = <)^e^-«=^, (42) 

where R^^^ is a normalization constant and Xbd is extracted from the numerical solutions in the regime where it is 
approximately constant. These extracted values for the solutions with two kernels as well as various values of as are 
found in Table |TT1 The table also shows that with changing as the exponent is relatively constant for the LO kernel, 
once again nonlinearities in the exponent appear for the Bessel kernel. Reported exponents are averaged from values 
of dipole size r = 10~^ 10^ because Abd does vary slightly with dipole size. 



VI. INCLUDING THE RUNNING COUPLING 

Turning now to the running of the QCD coupling, ctg is taken to be as{r) = ^^^^^ i — j where b = —j^^, "-/is the 
number of active flavors and A = 0.246 GeV is used. In the infra-red regime, the coupling was frozen when r > Tcut 
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FIG. 13: Graphs of the natural logarithm of the black disc cross section versus rapidity at various fixed dipole sizes. The 
dashed line represents the solutions obtained with Bessel kernel while the solid line represents the solutions with the LO kernel. 



which is defined as as{rcut) = 0.3. As it is well known the BK equation without the impact parameter is not very 
sensitive to the way the coupling is regularized. This is because the amplitude is saturated for all the large values of 
the dipole size from the inverse of the saturation scale to infinity. In the case with impact parameter however, there 
are contributions from the large dipole regime which spoil this self-regularizing behavior. In this case there is a large 
sensitivity to the regularization scenario for the running coupling. 

There are two different schemes for including the running coupling in the BK equation [65| . [6^ . In addition to 
these two scenarios we will use also the so-called parent dipole scheme, where the coupling depends on the size of the 
external dipole, that is xqi. This scheme is convenient to use with the Bessel function kernel. We have also evaluated 
the solutions using the prescription proposed in [65j 



iVcas(a^oi) 
27r2 



"-01 



1 



'-02-^12 



«^(^02) 

^02 \Q^s(a^i2) 



1 



n2 



Q^s(a^l2) 



02) 



(43) 



Since it is not clear at the moment how to use this scheme with the Bessel function kernel we will use it only with the 
LO kernel. The scheme dependence between two prescriptions 65], 66] originates from the choice of the subtraction 
point. The scheme by [1^ was shown to agree with the scheme [65j by the calculation of the appropriate subtraction 
corrections. In this paper we have not evaluated the scheme (66l |. as we have found that in order to achieve the desired 
accuracy for the solution with impact parameter within this scheme takes considerably longer time. 

We first shall show the results with the running coupling without the impact parameter. The running of the coupling 
has the effect of slowing down the evolution of the scattering amplitude as seen in FigfMl The difference between the 
LO and the modified kernel with running coupling is rather small. This can also be seen in Fig IT5] which shows the 
saturation scale of the two kernels with running coupling which are extremely close to each other. The dependence 
on the saturation scale with respect to the rapidity is as in [s^l 



A exp 



2x(7c 



■Y 



X (7c) 
2&7cX(7c) 



1/3 



Ye 



(44) 



where 



-2.338. Here the second term involving is numerically non-negligible for the rapidities we consider. 



In terms of numbers the coefficients above give — X^e^'^^^ 5.4Fs have found that the LO saturation scale 

1 1 

with running coupling and the parent dipole size prescription is = g'^-*^' which is very similar to the one 
given by the analytical value. The running coupling with prescription (j43|) has also been run and found to have a fit 



of 



which is closer to the value given by (jH]). 



In the scenario with impact parameter we find quite different behavior of the solution. As is seen in Fig. 16(a) the 



evolution of the running coupling (with parent dipole scheme) is actually very fast in the small dipole region, and 
it is much faster in the large dipole region. This is obvious since in the large dipole region the coupling is fixed at 
tts — 0.3 which yields approximately three times as fast an evolution versus the case where IXs = 0.1 is fixed. It can 
be seen there are box effects beginning to manifest in the running coupling case due to the frozen coupling evolving 
very quickly in the large dipole regime and reaching the box. 

It can be seen in Fig.[T7]that the dependence of the saturation scales on the rapidity is now again almost exponential. 
In this case we can extract the exponents by fitting exponential forms in the rapidity as we did for the fixed coupling 
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(a)Solid line: LO kernel with fixed as = 0.1; dashed line: LO (b)Solid line: LO kernel with = as{r); dashed lino: modified 
kernel with Us = as{r) kernel Us = as{r) 

FIG. 14: Graphs of the scattering amplitude versus dipole size for the case without impact parameter. The left graph illustrates 
how much slower propagation is due to the running of the coupling and the right graph shows the modified and LO kernels 
compared to each other in the case of the running coupling. Each line corresponds to the rapidity increasing in rapidity in 
intervals of AF = 10 up to F = 100. 
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FIG. 15: Graph of the saturation scale of the LO kernel (solid line) and the modified kernel (dashed line) with running coupling 
and no impact parameter dependence. 



case (table Hill) . Note that the definition of the exponents are now different than in the previous section. Here, we 
took Qs ~ exp{XsY),QsL ~ exp(— A^y), i?s ~ exp(ABy). The reason that the dependencies are almost exponential 
is due to the large sensitivity to the infrared and the fact that the coupling is frozen. In that case the solutions behave 
as almost with the fixed running determined by the freezing value. 

Similar pattern is found in the case of the running coupling with the scenario (|43p . The only difference is in the 
small dipole regime where the evolution is slightly slower than that of the parent dipole scheme. This can be seen by 
comparing the extracted exponents in Table Hill 
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FIG. 16: Left (a) dipole scattering amplitude as a function of the dipole size for fixed impact parameter and angle. Solid lines: 
fixed coupling with cts = 0.1, dashed lines: running coupling with the parent dipole scheme. Right (b): saturation scale as 
function of rapidity for the LO kernel the fixed coupling =0.1 (solid line) and running coupling (dashed line) case. 
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(a)Saturation scale at small dipole size (b)Saturation scale at large dipole size (c)Black disc radius in impact parameter 

FIG. 17: Small dipole saturation scale, large dipole saturation scale and the black disc radius for the case of the running 
coupling within the parent dipole scheme. The solid lines are for the LO kernel and dashed lines are for Bessel kernel. 



The behavior observed is of course something that has been analyzed before, in the context of the linear BFKL 
with running coupling ^6^. In particular it was observed that, the BFKL solution shows the tunneling scenario, where 
at some value of rapidity the solution is completely dominated by the infrared region. Strictly speaking we are not 
observing the tunneling scenario here, due to the fact that we have chosen our initial conditions to be concentrated 
around rather large dipole sizes where the coupling is already large. Rather, our solutions are completely dominated 
by the large coupling values and hence the saturation scale has nearly exponential dependence on rapidity. It will be 
interesting to analyze the solution for the initial conditions which are located in the small dipole regime to see if the 
tunneling occurs here. 

We have also evaluated the dipole cross section coming from the black disc regime in the running coupling scenario. 







Al 


As 


AsD 


LO Kernel q^'™' 


0.30 


1.68 


0.60 


0.65 


LO Kernel q,'-^"'' 


0.29 


1.68 


0.64 


0.68 


Bessel Kernel ai^'^^ 


0.22 


1.42 


0.24 


0.32 



TABLE III: Summary of extracted evolution exponents with impact parameter for running coupling case. PD means parent 
dipole prescription, Bal means prescription ((43}. 
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(a)Black disc cross section for the fixed coupling Os =0.1 
case (solid line) versus the running coupling case (dashed 
line), both for the LO kernel 
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2 4 6 8 10 12 
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(b)Black disc cross section for case of running coupling 
with the LO kernel (solid line) versus the running 
coupling case with the Bessel kernel(dashed line) 



FIG. 18: Dipole cross section. Contribution from the 'black disc'. 



and we parametrize it in the form 

aBD = 27TRlj,ix, Y) = 2TTR%e^^^^ . (45) 

The extracted value for the exponent \bd is shown also in Tabl dllll Again the black disc cross section increases 
very fast due to the large value of the coupling in the region of freezing. We have also compared the solutions in the 
case of the LO and Bessel kernel, the results are shown in Fig. [181 Since the coupling is relatively large, the differences 
between the evolution with LO and Bessel kernels are more amplified. 



VII. CONCLUSIONS 



Let us summarize the most important points of our analysis. 

In the case of the solution with the LO kernel, the extracted exponents of the saturation scales and the black disc 
radius are consistent with values obtained from the boundary method. The peaks in impact parameter and in dipole 
size can be very easily understood from arguments based on the conformal symmetry. The relative strength of the 
evolution of different saturation scales follows as well from the arguments on the conformal symmetry. In particular 
the black disc radius has an expansion rate which is twice slower than that of the saturation scale for small dipoles. 

For the running coupling scenario, in the case of the solutions with impact parameter we no longer observe the 
self-regularizing behavior of the nonlinear equation. This is of course due to the increased sensitivity to the large 
values of the dipole size. Rather, for the initial conditions chosen one observes strong dependence on the details of 
the regularization, and basically the exponents of both the saturation scales are dominated by the largest value of the 
coupling. This could be tested in more detail by choosing different initial condition, nevertheless one can expect that 
for the sufficiently large rapidity, the solution becomes regularization-sensitive, much like it was observed in earlier 
simulations. 

The cuts on the large dipole sizes, introduced in the form of the modified kernel have in general very small effect 
for the case of the impact parameter independent kernel. For the case with the impact parameter they are no longer 
negligible and reduce the exponent by about 25% for coupling of = 0.1. It is important to note that, the modified 
kernel we have chosen does not account for all the type of kinematical cuts, and therefore other cuts, on the small 
dipole sizes should be included similarly to what was done in fsoj. One could expect therefore even stronger effect in 
this case. 

We therefore conclude that the observed self-regularizing behavior of the local BK equation with the running 
coupling and almost complete insensitivity to the other NLO corrections appear due to the simplified assumption 
about the impact parameter independence. 

In a broader perspective, it will be interesting to perform the analysis with full NLO kernel or the more correct 
form of the kinematical cuts, as well as introduce effectively confinement effects. It is also vital to analyze the impact 
of the corrections which go beyond the mean field approximation [62,, .68l - l70l | . 
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